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Abstract 

We find a (nonlocal) gauge where the wavefunction renormalization con- 
stant does not get any corrections for all momenta in the hard-dense loop 
approximation. In this gauge, we solve the Schwinger-Dyson equations for 
the diquark condensate in dense QCD to calculate the Cooper pair gap. We 
determine not only the exponent but also the prefactor of the gap in a gauge 
independent way. We find that the higher order corrections increase the gap 
only by about 1.6 times to the leading order gap at Coulomb gauge. 

PACS numbers: 12.38.Aw, 12.38.Mh, 11.15.Ex 
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Typeset using REVTgX 



Matter at extreme density is known to be a color superconductor [1], quarks forming 
Cooper pairs to open a gap at the Fermi momenta. The properties of color superconductor 
are intensively studied recent years at asymptotic density [2], where quarks interact weakly 
due to asymptotic freedom in quantum chromodynamics (QCD). However, the study beyond 
the leading order has been hindered by calculational difficulties. 

Since superconductors are characterized by the Cooper-pair gap or the minimum energy 
to excite a pair of a quark and a hole at the Fermi momenta, it is quite important to 
determine the size of the gap accurately. By solving the Schwinger- Dyson (SD) equations, 
the gap at the asymptotic density has been found to be 

Ao^expf--) (1) 

9 s \ 9s J 

where the l/g s power in the exponent is due to the long-range interaction in dense QCD. 
The constant c in the exponent is gauge-independent and is determined by the long-range 
interaction mediated by the magnetic gluons, unscreened at high baryon density. In the 
hard-dense loop (HDL) approximation, it is known to be c = 3ir 2 /\/2 [3]. The short- 
range interaction mediated by the Debye-screened electric gluons contributes only to the 
prefactor. Among the negative fifth power of the coupling in the prefactor, the third is due 
to the screened electric gluons, while the magnetic gluons contribute the remaining. In the 
leading order analysis of the SD equation [4,5] the numerical prefactor b is found to be 




2 V * e 3,/ 2+ i_ (2) 



where Nf is the number of quark flavors and £ is the gauge-fixing parameter. 

The result Eq. (2) is gauge-dependent, since all subleading contributions that lead to 
logarithmic divergences as in the BCS superconductivity are not taken into account. Being 
a physical observable, the gap A should be gauge-independent. However, if one sums the 
contributions partially, the result is often gauge-dependent. We therefore need to calculate 
the subleading corrections to obtain the gauege-independent gap. Some of the subleading 
corrections like the finiteness of the quasi-quark life in medium [6], the running effect of the 
strong coupling [7], and the quark self-energy [8,9] have been studied. 

In this paper we determine accurately the contributions to the prefactor, coming from 
the vertex corrections and the wavefunction renormalization for quarks, since these are the 
remaining gauge-dependent contributions to the prefactor and, if summed up, we should get 
a gauge-independent prefactor. 

We first try to find a (nonlocal) gauge [10] where the quark wavefunction is not renor- 
malized for all momenta, Z(p) — 1. In this gauge we then solve the SD equations to find 
the Cooper pair gap, paying special attention to the Ward-Takahashi identity. This way of 
calculating higher order corrections in the SD analysis has been proved extremely useful in 
dynamical mass generation of (2+l)-dimensional quantum electrodynamics and in others, 
since the higher order corrections vanish exactly [11,12]. 

Unlike ordinary electron superconductors, without additional interaction the quark-quark 
scattering is attractive in the color antitriplet channel, where the color flux energy is low- 
ered. By the Cooper theorem, diquark operators of opposite momenta therefore develop a 
condensate in quark matter, breaking the color gauge symmetry, 
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x)ip c (x)) = K(p F ), 



(3) 



where the charge-conjugated field is defined as (ip c )i (x) = Cij"4>j(x)- The matrix C satisfies 
C _1 7 M C = — 7J and i, j are Dirac indices. 

To calculate the condensate, we introduce a Nambu-Gorkov field ^(x) = (^(x),^ c (x)) T . 
The inverse propagator for the Nambu-Gorkov field is then given as 



a(p) [(Po + /-Ot + Kp) f>\ -A(p) 

- 7 °At(p) 7 o a(p) [(po _ ^ )7 o + 6(p) #] 



(4) 



where a(p),b(p) denote the wave function renormalization constants and A(p) denotes the 
Cooper-pair gap at a momentum p. Notice that the time and spatial components of quark 
wavefunction renormalize differently since the Lorentz symmetry is broken in quark matter. 

We study the Schwinger-Dyson equation for the Nambu-Gorgov propagator, following 
the notations in [4]. The gap equation in dense QCD takes a following form (See Fig. 1); 
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The first term in Eq. (5) comes from the magnetic gluons and gives the leading contri- 
bution. The constant term, In b, is due to the Debye screened electric gluons and the gauge 
fixing terms. The higher-order terms, denoted as i] and (, are in general given in powers 

1 [ A n ~\ m [ 

of coupling constant and energy as g s — ^ In 



_ i n our calculation q ~ A. The 



higher-order corrections are in general suppressed by powers of coupling constant. However, 
since the gap depends on the coupling, the logarithmic corrections might make the higher- 
order corrections comparable to the leading term. We examine the two-loop corrections to 
the gap equation in detail. 

In general, the SD equations are infinitely coupled equations for the 1PI functions. How- 
ever, some of them are related by Ward-Takahashi (WT) identities. The identity we are 
going to use in solving the gap equation is one that relates the quark two-point functions 
with the vertex functions, 



(6) 



= (d%(z)^(x)^(y)) - 5{z - x) (TV(s$(i/)) + 5(z - y) (^{x)^y)T a 
The identity reads in the momentum space as 



T a o(p) [(po + /i)7° + Hp) #] - <P') [(Po + 0)7° + Hp') / 
= (p-p')^(p,p')T a + T a (p,p';-p-p'), 



rpa 



where 



T a (p,p';k)5(k + p + p') = I , 

Jz,x,y 



i(z-k+x-p+y-p') 
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(7) 



(8) 



The one-loop vertex correction has two parts (See Fig. 2). 

As in quantum electrodynamics, the first part (Fig. 2a) is related to the correction to 
the wavefunction renormalization constant as 
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(p - p'),K^{p,p') = a(p) f(po + /i) 7 ° + 6(p) f] ~ a(p') \(p' + /i) 7 ° + btf) f 



where we suppressed the color indices. We see that A*"^ = 7^ for the nonlocal gauge where 
the wavefunction renormalization constants are a(p) = 1 = b(p). (This result is similar 
to the condition obtained by Kugo and Mitchard in the SD analysis of chiral symmetry 
breaking in QCD, using vector Ward identities [13].) Therefore, if we find a gauge where 
a(p) = 1 = b(p), not only the corrections from the wavefunction renormalization but also 
those from the vertex corrections from Fig. 2a are absent. This simplifies significantly the 
calculation in the higher order corrections to the Cooper-pair gap, since the vertex correction 
to the gap, coming from the second diagram, Fig. 2b, does not contribute to the prefactor, b, 
in Eq. (1). To see this, we note the leading contributions to the vertex comes from when both 
the internal and the external gluon lines are magnetic. Furthermore, the main contribution 
to the gap occurs when the external lines carry small momenta. For small external gluon 
and quark momenta (q, \p\ — /i, — /i ~ A), the vertex correction due to the diagram 
Fig. 2b is given as, using the high density effective theory of QCD [5], 
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where c\ and C5 are 1 + (A 2 //^) while all other q's are O (A 2 //^). The screening mass M 
is given as g s j2^jNf/(27r) for Nf light quarks in the hard-dense-loop (HDL) approximation. 
Therefore, we find upon integration the correction becomes 



^~^ s 3 7o^T a (^ 2 i n( A/^), 



(10) 



which is indeed negligible compared to the constant term, In b, in the gap equation Eq. (5). 

Now, we look for a nonlocal gauge where the wave function constants remain unrenor- 
malized, a(p) = 1 = b(p). At this gauge, the vertex is bare, A M = 7^, and the SD equation 
leads to an equation for the (nonlocal) gauge-fixing parameter, given as 



d\D, u {q - p)Tr ( 7 VA+ 7 VA; 



(gb + |g1 ~ V) 
9o 2 -(|g1-^) 2 -A 2 



0, 



(11) 



where the antiquarks are projected out by the positive and negative energy state projectors, 
defined as 



(9) 




(12) 



The (Higgsed) gluon propagator is given in the HDL approximation as 1 



For two flavor color superconductors, some of the gluons are not Higgsed. However, since the 
Meissner effect is subleading, its inclusion does not affect our calculation. 
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D, l/ (k) = AO$ + BO$ + CO$, 
where in the weak coupling limit, \k \ <C \k\, 

k 



(13) 



A 



B 



C 



(14) 



k" + M 2 A + irM 2 \k 4 \/2 k\ + k 2 + 2M 2 ' kj + k 2 ' 

where M is Higgs-like gluon mass ~ g a \xj (2%) [4]. The polarization tensors are defined as [4] 

(D = p x + (u-kf 



_pu 

(u ■ k) 2 — k 2 ' 

(2) = p±_ (l) ) (3) = p|| ) 



(15) 



where = (1,0) and 



P x = 

kp}z v 


kjj,u u + m^A; 
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(16) 



C^Tr ( 7 7 ^A+ 7 7 ^A;) = -2(1 + t) 
0%Tr ( 7 V A+ 7 ° 7 ^A;) = 1 - t 
O^Tr( 7 VA+7VA;) = -(l-t) 



q 2 +p 2 - qp(l + t) 

g2 _|_ ^2 _ 
(P + ?) 2 



g2 _|_ ^2 _ 2pg£ 



(17) 



where p = g = |g|, and t — p ■ q/(pq). 

Now, we assume the (nonlocal) gauge fixing parameter £(k) has only temporal depen- 
dence, that is, ~ £(k ), and would like to have (11) satisfied already at angular inte- 
gration. This gives 



D + E + F = 



(18) 



where 



D = q I ^dtA(p 4 - q 4 ,p- q) 



q2 _|_ p2 _ iqpt 
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(2/i) ; 



(q - p) 3 + M 2 A + nM 2 \p 4 - q 4 \/A 
E = q 2 J dtB(p 4 -q 4 ,p-q)(l-t) « 1 

F = q 2 J ^dtC(p 4 - q 4 ,p- q) 
(2/i) 2 



g 2 + p 2 — 2pgt 



-£ln 



(p - g) 2 + |p4 - 94 1 



(19) 
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In this calculation, we take p,q ~ since the quarks are near the Fermi surface. We find 
the solution of Eq (18) 



2l n (2M) 3 



3 M^A+7rM 2 |p4-94|/2 



1 

3' 




(20) 



with |p4 — g 4 1 ~ A. 

Plugging Eq. (20) into the leading-order expression for the gap Eq. (2), we get 




(21) 




In conclusion, we find that the wave-function renormalization vanish at a nonlocal gauge- 
fixing parameter £ ~ 1/3. By the Ward-Takahashi identity, the QED-like vertex correction 
vanishes at this gauge as well. Since the remaining one-loop vertex corrections are suppressed 
for small external momenta, which has a major contribution to the gap equation, both the 
wave-function correction and the vertex correction to the prefactor of the gap vanishes at 
this nonlocal gauge. As the gap is independent of the choice of gauge, we conclude that the 
next-to-leading order corrections due to the vertex renormalization and the wavefunction 
renormalization increase the leading-order gap at the Coulomb gauge by about two thirds. 
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FIGURES 



■ A ■ 

FIG. 1. The triangle denotes the Cooper-pair gap and the blobs higher-order corrections. 




FIG. 2. The solid line denotes quarks and the curly lines gluons. 
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